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ABSTRACT 

Astrometry can bring powerful constraints to bear on a variety of scientific questions about neutron 
stars, including their origins, astrophysics, evolution, and environments. Using phase-referenced ob- 
servations at the VLBA, in conjunction with pulsar gating and in-beam calibration, we have measured 
the parallaxes and proper motions for 14 pulsars. The smallest measured parallax in our sample is 
0.13 ± 0.02 mas for PSR B1541+09, which has a most probable distance of 7.2+};^ kpc. We detail 
our methods, including initial VLA surveys to select candidates and find in-beam calibrators, VLBA 
phase-referencing, pulsar gating, calibration, and data reduction. The use of the bootstrap method to 
estimate astrometric uncertainties in the presence of unmodeled systematic errors is also described. 
Based on our new model-independent estimates for distance and transverse velocity, we investigate 
the kinematics and birth sites of the pulsars and revisit models of the Galactic electron density distri- 
bution. We find that young pulsars are moving away from the Galactic plane, as expected, and that 
age estimates from kinematics and pulsar spindown are generally in agreement, with certain notable 
exceptions. Given its present trajectory, the pulsar B2045— 16 was plausibly born in the open cluster 
NGG 6604. For several high-latitude pulsars, the NE2001 electron density model underestimates the 
parallax distances by a factor of two, while in others the estimates agree with or are larger than the 
parallax distances, suggesting that the interstellar medium is irregular on relevant length scales. The 
VLBA astrometric results for the recycled pulsar J1713-I-0747 are consistent with two independent 
estimates from pulse timing, enabling a consistency check between the different reference frames. 
Subject headings: astrometry — pulsars: individual (B0031-07, B0136-I-57, B0450-18, B0450-I-55, 
J0538+2817, B0818-13, B1508+55, B1541+09, J1713-f'0747, B1933+16, 
B2045-16, B2053-H36, B2154+40, B2310+42) — stars: distances — stars: kine- 
matics — stars: neutron 



1. BACKGROUND AND GOALS 

Neutron stars are exotic laboratories for some of the 
most extreme physics in the Universe. Precise and accu- 
rate measurements of the position, proper motion, and 
parallax of neutron stars (NS) can be obtained by as- 
trometry, and such measurements can be exploited to 
address a variety of scientific questions. 

For example, measuring the position of an object in two 
different coordinate frames permits the very fundamen- 
tal operation of tying the two reference frames together. 
Since they are compact sources and are accessible to as- 
trometry with different techniques and at different wave- 
lengths, NS are particularly well suited to such reference 
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frame ties. Specifically, the optical counterpart of the ra- 
dio pulsar J0437— 4715 may provide a frame-tie between 
the optical reference frame and the radio- defined Interna- 
tional Celestial Reference Frame (ICRF; iMa et aLlfTOOSh 
for the Space Interferometry Mission. 

Precise proper motions for NS allow them to be traced 
back to their birth sites in massive stellar clust ers, and to 
associations with runaway stars jH oogcrwcrf et al.|[2000l : 
IVlemmings et al.l l2004': Cha tteriee et al.1 12005'). For very 
young objects, associations with their progenitor su- 
pernova remnants can be verified or refuted, leading 
to independent age estimates fo r both the N S and the 
supernova remnant itself (e.g.. iGaensler fc Fr ail 200d 
Thorsctt ct al. 2002: Migliazzo ct al." 20021 iKramer et all 
2003; Blazek et al. 2006: Zeiger et al. 20081). 

Combined with estimates for their distances, the 
proper motions of NS also lead to velocity estimates. The 
high velocity tail of the distribution implies that large 
kicks are imparted to proto-neutron stars during core 
collapse, and a number of mechanisms have been pro- 
posed for these iiatal kicks (e.g., [B urrows & Haves 1991; 
Janka fc MuelieilfTOQ^ IA rras fc Laiiil999. : Socrates et aLl 
20051) . The proper motions of magnetars (Hclfand et al.l 
; iDe Luca et "alll2008l : iKaplan et all [2008b ) can po- 
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tentially test the role of strong magnetic fields in such 
kicks, while precise comparisons of the proper motion 
and spin axis vectors of pulsars (iJohnston et al.l l2005t 
iNg fc Romallill2007t iKaplan et al.ll2008af l constrain the 
degree of rotational averaging in hydrodynamic kicks 
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(e.g. lLai et alll2001[ ). 

When a parallax measurement is possible, a model in- 
dependent estimate is obtained for the distance and veloc- 
ity of the NS. Each such measurement calibrates global 
models of the G alactic electron density ([Taylor &: Corded 
119931 ICordes fc Lazio 2002, ,) , thus improving distance es- 
timates from pulse dispersion measure for the rest of the 
radio pulsar population, as well as probing the distribu- 
tion of electron density in the local interstellar m edium 
(;e.g.. lToscano et al.|[l99l IChatteriee et al.ll2001h . 

Observed thermal radiation from the NS surface can 
be used, in combination with a precise distance, to 
constrain the size of the photosphere, the NS ra- 
dius, and thus the Equation of State of matter at ex- 
treme pressures and densities fYako vlev fc Pethickll200 j : 
iLattimer fc Prakash,2004) . For radio pulsars, uncertain- 
ties in the magnetospheric emission restrict such an exer- 
cise to very young an d hot objects (e.g., PSR B0656-I-14, 
iBrisken et al] l2003bf) . while isolated NS which are X- 
ray bright and radio quiet pose a chall enge for optical 
astrometry (e.g., RX J18 56. 5-3754, Kaplan et all 120021 : 
IWalter fc 'Lattimeill2002j ). 

Astrometry can thus bring powerful constraints to bear 
on a variety of scientific questions about neutron stars, 
including their origins, astrophysics, evolution, and envi- 
ronments. The primary obstacle is the difficulty of such 
astrometric observations. As in most other astrometric 
applications, multiple position measurements of a given 
target over some time span are the basic observables from 
which astrometric parameters are derived. For individ- 
ual objects, rather than an ensemble of stars over the 
entire sky, such measurements are necessarily relative in 
nature, but absolute positions can be inferred from mea- 
surements relative to sources that define the ICRF. Over 
time, repeated measurements of the position 9 allow a 
proper motion pi to be derived. For a precise proper mo- 
tion, the primary consideration is a long time baseline, 
limited by the stability of the reference frame and the 
variability or motion of the frame-defining sources. Fi- 
nally, with enough astrometric precision, a trigonometric 
parallax tt may also be measurable for NS. The primary 
consideration for such measurements is appropriate sam- 
pling over the Earth's orbital phase, not just a long time 
baseline. 

Neutron stars emit over a broad range of frequencies, 
and astrometric observations of NS have been conducted 
at bands froi n radio (~ 10^ Hz) t o X-rays (~ 10^'^ Hz). 
For example, I Kaplan et al.l (|2007t l have used optical ob- 
servations with the Hubble Space Telescope to determine 
a proper motion /i = 107.8 ± 1.2 mas yr~^ and a paral- 
lax TT = 2.8 ± 0.9 mas for RX J0720.4-3125. However, 
only a few NS can be observed at optical or infra-red fre- 
quencies. In X-rays, the resolution that can be achieved 
with the c urrent generation of te lescopes is a limiting 
factor, but I Winkler fc Petrd (j2003) have used the Chan- 
dra X-ray Observatory to measure a proper motion of 
165 ± 25 mas yr'^ for the RX J0822-4300, in the center 
of the Puppis A supernova remnant. 

The majority of known NS are radio pulsars, and 
pulse timing is routinely used to refine their positions 
and proper motions. A subset of recycled ("millisec- 
ond" ) pulsars have rotation rates that are stable enough 
to permit sub-milliarcsecond astrometry based on pulse 
time of arrival. As a particularly interesting example. 



IVerbiest et al] (|2008[ ) use pulse timing at Parkes to in- 
fer a precise distance for the binary pulsar J0437— 4715 
based on the observed orbital period derivative Pf, and 
the assumption that the intrinsic orbital period deriva- 
tive is dominated by the emission of gravitational ra- 
diation. The measurement leads to an estimate for 
the pulsar mass and an upper limit on the varia- 
tion of the gravitational constant G. However, most 
pulsars do not have such stable rotation, particularly 
when they are young, and Very Long Baseline inter- 
ferometry (VLBI) has usually been utilized to deter- 
mine their astrometr ic parameters. Su ch efforts have 
a long history (e.g., iGwinn et al.lll986l ). but have be- 
come much more feasible with the Very Long Base- 
line Array (VLB A), which provides full-time, dedicated 
VLBI capabilities with identical antennas, allowing good 
control of systematic errors and leading to many re- 
cent pa rallax measurements (e.g., iFb malont et al.' '199!^ 
Brisken et al. 2000; Chattcri ee"et al.l [2001: Briskcn et all 
2002t IChatteriee et al.l l2004l . 

However, not only are most pulsars faint, but the most 
interesting categories, namely the youngest pulsars and 
the recycled ones, appear to be disproportionately faint. 
Pulsar gating (which accumulates signal only during the 
predicted on-pulse periods for pulsars) is now routinely 
employed to boost the signal-to- noise ratios (S/N) for 
VLBA astrometric observations, using contemporaneous 
pulse timing to reliably predict the pulse phase at the cor- 
relator. To make matters worse, however, it is not just 
a lack of sensitivity but systematic errors contributed 
by the ionosphere and the troposphere which are the 
primary impediments to sub-milliarcsecond astrometry. 
Progress has been made, for example, with GPS-based 
ionospheric calibration schemes (Walker fc Chatterjee 
1999, VLBA Scientific Memo^" 23) and with wideband 
ionos pheric calibration tec hniques (Brisken et al.l l2000l . 
|2002|) . In-beam calibration (jFomalont et al.lll999| ). using 
a faint source in the same primary be am as the target 
source, has also proved to be effective (jChatteriee et al.l 
120011 l2005h . 

Here we present a set of results from a large astrometric 
program to determine the proper motions and parallaxes 
of radio pulsars. The project involved a large scale initial 
survey at the Very Large Array (VLA) to identify target 
pulsars and detect in-beam calibrators, as discussed in 
|j2l followed by multi-epoch observations with the VLBA 
(SjS)). We provide details of the calibration and data re- 
duction in m and our use of the bootstrap method to 
estimate uncertainties is discussed in [J5l as a template for 
ongoing and future VLBA astrometry programs. Some 
implications of our results are explored in SJH specifically 
for NS birth sites. Galactic electron density models, and 
reference frame ties. Finally, our conclusions are sum- 
marized in fjTl 

2. PULSAR SELECTION AND CALIBRATOR 
IDENTIFICATION 

Astrometric measurements with the VLBA are con- 
ducted in the International Celestial Reference Frame 
(ICRF, iMa et al.lfl998l : iFev et al.l[200l . Each pulsar of 
interest thus requires the identification of suitable nearby 
calibrators that tie the observations to the ICRF. 

See |http : //www . vlba ■ nr ao . edu/memos/sci/| 
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In the northern hemisphere, the density of ICRF 
sources and secondary calibration sourc es with suffi- 
cientl y accurate posi tions in the ICRF (iBeaslev et alj 
I2002t iFomalont et al.l 120031 : iPetrov et al.ll2005[ ) is about 
0.1 deg~^, implying typical calibrator-target separations 
of '^2°, although the nearest calibrators can be as much 
as 5° distant in some cases. Astrometric errors due to dif- 
ferential propagation effects through Earth's ionosphere 
and troposphere scale with calibrator separation (e.g., 
IChatteriee et al.l 12004V requiring smaller angular sepa- 
rations for higher astrometric accuracy. The detection 
and use of compact sources considerably closer to the tar- 
get than those in the existing catalogs offers a straight- 
forward way to improve astrometric precision. Unfor- 
tunately, only a small fraction of the many sufficiently 
bright point-like sources typically detected by the VLA, 
even in its most extended (highest resolution) configura- 
tion, are compact enough (< 10 mas) to serve as VLBA 
phase calibrators, owing either to intrinsic source size 
or to scattering by the intervening interstellar medium 
(ISM). 

In our previous work, the search for such nearby com- 
pact calibrato rs has led to the develg pr nent of "in-beam 
calibration" (jFomalont et al.l Il999t IChatteriee et al.l 
I200I . Here a compact source is identified within the 
same field of view (i.e., primary beam) as the target 
source. For the 25-m dishes of the VLA and VLBA, 
this implies an angular separation less than ~25' between 
the target and the in-beam source for observations at the 
upper end of our band (~ 1.7 GHz). Such sources are 
likely to be weak, requiring the use of a stronger, more 
distant phase referencing source for initial calibration at 
each epoch, but the position of the target is thereafter 
referenced to the in-beam calibrator position. In-beam 
calibration has the dual benefit of allowing the calibra- 
tor and target to be observed simultaneously, thereby in- 
creasing the sensitivity for both objects, as well as reduc- 
ing the angular extrapolation of the calibration solutions 
and eliminating the need to interpolate the calibration 
in time. 

2.1. Initial Pulsar Field Survey 

An initial sample of 63 pulsars was chosen from the 
roughly one thousand pulsars known at the inception of 
this project in 2001. Pulsars were selected based primar- 
ily on flux density, declination, and dispersion measure 
distance estimate, so that precision astrometry was plau- 
sible using techniques refined in our previous work. A few 
pulsars were included since they were of particular scien- 
tific interest (e.g., the recycled pulsar J1713-I-0747) even 
though they were expected to pose technical challenges; 
sample completeness was not a consideration. 

The field of view around each of the 63 selected pul- 
sars was observed with the VLA between 2002 Febru- 
ary 7 and 9 (project code AC629). Given the similar 
25-m diameters of the VLA and VLBA dishes, they have 
well matched primary beams at the same observation 
frequencies. Each field was observed for ^^15 minutes at 
1.4 GHz with the VLA in its highest angular resolution 
A-configuration (~1'.'4 at 1.4 GHz). A spectral-line ob- 
serving mode was used, with a 23 MHz spectral window 
or bandpass, and the correlator output was dumped ev- 
ery 5 s (as opposed to the more standard 10 s). The use 
of a spectral-line mode and short correlator dump times 



were both motivated by the desire to image the full field 
of view with minimal distortions due to time- average and 
bandwidth smearing (|Bridle fc Schwablfl999l) . 

In order to detect potential in-beam calibrators, as 
well as bright confusing sources within the primary 
beam and beyond, a low-resolution image with a diam- 
eter of 1?1 was searched for sources. Small facets sur- 
rounding each source were then imaged and deconvolved 
at the full resolution, employing polyhedral imaging 
(|Cornwell &: Perlevlll992f ). For each field, self-calibration 
was used to improve the image dynamic range, and if a 
substantial improvement was obtained in the dynamic 
range, the source-finding process was repeated in an ef- 
fort to identify any weaker sources that might have been 
missed in the first iteration. 

Each detected source was inspected for compactness 
and cataloged with its brightness, position, and angu- 
lar separation from the pulsar. A total of 1060 sources 
stronger than 1 mJy were detected in the 63 fields; of 
these, 269 were within 25' of the target pulsar and were 
deemed compact enough (unresolved at 1'.'4 resolution) 
to warrant further investigation. At least one such candi- 
date in-beam calibrator was found for 62 of the 63 target 
pulsars, and an example is shown in Figure [1] In ad- 
dition, the position of each pulsar was measured to a 
precision of between 0'.'2 and 0'.'05, sufficient to enable 
the later VLBA observations. 

The angular resolution of these preliminary VLA ob- 
servations was still a factor of 100 below that of the 
planned VLBA observations. A second set of VLA ob- 
servations at 8.4 GHz was used to determine the sub- 
arcsecond structure and verify the compactness of the 
candidate in-beam calibrators, as well as to determine 
their spectral indices. Each of the 269 candidate in-beam 
calibrators was observed in individual pointings on 2002 
April 6 with the VLA A-configuration at 8.4 GHz, which 
provides a resolution '--^0'.'2. Each source was observed for 
approximately 3 minutes in continuum mode with two 
50 MHz spectral windows. Sources that were resolved 
were rejected, leaving 97 compact sources as potential 
in-beam calibrators with the VLBA, with at least one 
for all but eight of the original 63 pulsars. The refine- 
ment process is illustrated in Figure [TJ and details of all 
the candidate sources are archived online^^. 

2.2. Final Target Selection 

From the original sample of 63 pulsars, 34 were se- 
lected for observations with the VLBA, most of them 
with good candidate in-beam calibrators. A few pulsars 
were judged to be strong enough and to have a cata- 
loged VLBA calibrator that was close enough that phase 
referencing, in combination wi th wideband ionospheric 
calibration (jBrisken et al.ll200^ . would provide accurate 
astrometry. These pulsars will be discussed in future 
work. 

In the first epoch of VLBA observations, most of the 
selected pulsars were found to have at least one suit- 
able in-beam calibrator, as illustrated in Figure [TJ How- 
ever, the in-beam candidates for 7 pulsars were either 
not detected or were too heavily resolved to be useful. 
In addition, 3 pulsars were not detected in first epoch 
VLBA observations due to scattering or poor calibration 

|http : //www ■ astro ■ Cornell . edu/~shainl/psrvlb/ 1 
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Fig. 1. — An illustration of the selection and refinement process for in-beam calibrator candidates for PSR B2045— 16. Note the change 
in field of view for each row. Top row: Images at 1.4 GHz, VLA A-configuration, showing the pulsar and two compact in-beam sources. 
The beam size is ~1'.'4, and the image RMS noise level is ~0.5 mjy/beam. Besides the pulsar, 14 sources were detected in the field, and 
6 were found to be relatively compact. Middle row: Follow-up observations at 8.4 GIfz, VLA A-configuration. The beam size is ~0'.'2, 
and the image RMS noise level is ~0.3 mJy/beam. The pulsar was not observed at this high frequency, but all 6 targets were; 2 of the 6 
(shown here) are observed to be compact. Bottom row: First epoch VLBA observations at 1.4 GHz. The beam size is 14 X 5 mas, and 
the ungated image RMS noise level is ~0.3 mJy/beam. The pulsar and both in-beam calibrator candidates are detected. Note that for 
plotting purposes only, the flux densities of the in-beam candidates at the VLA, 1.4 GHz are not corrected for primary beam attenuation, 
and so apparently increase in the pointed VLBA 1.4 GHz observations, even though in practice some of the source flux density is probably 
resolved out. 

at lower declinations; these were discarded from our ob- 
serving list, and substituted with other targets. 

The process outlined here was optimized to identify 
pulsars where a parallax measurement seemed plausible 
based on our previous experience. No attempt was made 
to preserve sample completeness or to image each pos- 
sible in-beam calibrator candidate exhaustively, and as 
such, it is difficult to infer detailed statistical properties 
for the ensemble of faint, compact background sources. 
Of the 97 sources that were found to be compact in 
both 1.4 GHz and 8.4 GHz VLA A-configuration obser- 
vations, the vast majority (75 sources) had a negative 
spectral index {S^ oc z^"; a < 0), as expected for non- 
thermal sources. (The spectral index was determined 
after correcting measured flux densities at 1.4 GHz for 
primary beam attenuation.) Specifically, a < —0.5 for 
43 sources, and a < —1.0 for 16 sources. We empha- 
size that these results do not necessarily contradict the 



notion that com pact sources teii d to have flat spectral 
indices (e.g.. iPreston et al.lll983f ). Our sample is faint, 
and likely to include a mixture of different kinds of ob- 
jects, including some core-jet sources. Additionally, our 
two-point spectral index might straddle a peak in the 
spectrum, producing an apparent distribution of spec- 
tral mdices_tlmMshax^^ interpret. .Wall et al.i ()2005h 
and lFomalont et aTl ()20060 find broadly similar spectral 
index distributions for their source samples. 

3. VLBA OBSERVATIONS 

While pulsar astrometric observations have previously 
been conducted on a number of VLBI networks, the 
VLBA has several features that are particularly well 
suited to astrometry. The identical dishes of the VLBA 
simplify the calibration and ensure control of system- 
atic errors, and the identical primary beams (as well as 
the match with the primary beam of the VLA dishes) is 
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particularly advantageous for in-beam calibration. In a 
heterogeneous array, in-beam calibrators would have to 
be restricted to the field of view of the largest telescope 
in the network, which can severely restrict the chances 
of finding a suitable calibrator. Secondly, full-time oper- 
ation and flexible dynamic scheduling allow several short 
observations of each target through the year, timed to 
track parallax extrema and to avoid observations through 
the solar corona. Thirdly, the VLBA was designed for 
phase referenced observations with rapid slews and short 
dwell times, and provides repeatable phase stability, an 
essential component for astrometry. Finally, pulsar as- 
trometry in particular can take advantage of pulsar gat- 
ing, where the VLBA correlator effectively accumulates 
signal only when the radio pulse is predicted to arrive 
and discards the noise at other times. 

3.1. Observing Strategy and Scheduling 

Our choice of observing frequency balances compet- 
ing priorities. Radio pulsars typically have steep spectra 
(5*1/ oc a ~ —1.6), and thus are brighter at lower fre- 
quencies. Additionally, the larger field of view at lower 
frequencies enhances the likelihood of finding in-beam 
calibrators, and tropospheric phase fluctuations are re- 
duced at lower frequencies. On the other hand, the syn- 
thesized beam is broader, resulting in lower astromet- 
ric precision when sensitivity is the limiting factor (the 
case for many of the targets and calibrators here), and 
the phase distortions due to the differential ionosphere 
between the target and the calibrator increase rapidly 
at lower frequencies. An observing frequency around 
1.5 GHz was chosen as a compromise between these 
effects, as dem onstrated in our var ious p ilot projects 
(|Brisken et al.ll2000t[Chatteriee et al.ll2001[ ): the residual 
differential ionospheric effects remain the largest source 
of astrometric uncertainty in our program. 

In an effort to minimize systematic errors between 
epochs, all of the observations were set up in as identi- 
cal a manner as possible. Dual circular polarization was 
recorded simultaneously in four separate 8-MHz spectral 
windows. These bands were Nyquist sampled with 4-level 
quantization. We chose widely spread spe ctral windows 
to en able calibration of the ionosphere (|Brisken et al.l 
[2OOOI) for cases where in-beam calibration proved difficult 
and the pulsar was deemed bright enough to allow this 
method to succeed. The in-beam calibration technique 
described in the present work does not make use of the 
wide-spread bands. The initial set of spectral windows 
was centered on 1414, 1489, 1594, 1684 MHz, choices that 
were informed by presence of known satellite-generated 
interference. After routinely encountering strong radio 
frequency interference (RFI), the second spectral win- 
dow was changed to a central frequency of 1508 MHz. 
These spectral windows proved to be relatively RFI-free 
for most of the observations, although strong narrow- 
band RFI was occasionally observed around 1684 MHz, 
mainly at the Mauna Kea site. The VLBA pulse calibra- 
tion system was turned off during these observations in 
order to reduce the total system temperature by about 
3%. 

At each epoch, an observation cycled between the 
field containing the pulsar and its in-beam calibrator(s) 
and a VLBI reference calibrator. The VLBI reference 
calibrators, which are listed in Table [1] were drawn 



from the standard VLBI catalogs (iBeaslev et al.l 120021 : 
iFomalont et~ani2003l : iPetrov et al.ll2OO50 and are the ul- 
timate means by which the observations are tied to the 
ICRF. The calibrator was chosen primarily for proximity 
to the target pulsar, though strength and compactness 
requirements drove some choices. The pointing direction 
for a target field was set to be between the pulsar and 
its in-beam calibrator, with the restriction that the offset 
from the pulsar be no more than 6'; for some pulsars, two 
in-beam calibrators could be identified and the pointing 
was set to lie between all three sources in the field. Each 
observation spanned a duration of 2 hours approximately 
centered on the transit time of the source at the center 
of the VLBA. We used interleaved scans of 90 seconds 
on the reference calibrator and 120 seconds on the target 
field, resulting in a total integration time of about 1 hour 
on the target field after accounting for slewing overheads. 

We aimed to observe each pulsar for a total of 8 epochs 
over the course of 2 years, with epochs spaced approx- 
imately 3 months apart. Observations for a given pul- 
sar were conducted at a similar LST at each epoch in 
order to provide nearly identical spatial frequency sam- 
pling ("U — V coverage"). The observations were typi- 
cally scheduled within 2 weeks of parallax extrema. Be- 
tween 2002 June 27 and 2005 March 25, our observing 
programs (BC120, BC123) accumulated a total of 278 
observation blocks of two hours each. 

A few observations failed and had to be discarded due 
to the proximity of the Sun to the target pulsar in the 
sky, a circumstance whose severity was not fully appre- 
ciated by us when the program began. Turbulence in 
the extended solar corona and the solar wind can cause 
phase fluctuations (jSpangler et al.|[2"002f) which degrade 
astrometry. If the target field was within ~10° of the 
Sun, the images of the calibrators were severely distorted, 
but the solar corona and wind can affect 1.5 GHz obser- 
vations of targets as far away as 45° from the Sun. In 
order to avoid severe phase errors, future astrometric ob- 
servations should not be scheduled within about 30° of 
the Sun at 1.5 GHz. 

3.2. Pulsar Timing and Gating 

A priori information about the rotational phase of a 
pulsar can be exploited to correlate the observations only 
when the pulsar is "on" and the pulse is visible, while 
discarding data from the remainder of each pulse period. 
The result is a S /N improvement scaling as w ^ P/w for 
a pulsar with pulse period P and pulse width u;, with 
typical improvement factor ~3. The ephemerides used 
for gating the pulse signal were obtained from timing 
observations with the 76-m Lovell telescope of the Jodrell 
Bank Observatory. Typically, observing spans of a few 
months and overlapping the VLBA observation epoch 
were used to refine the pulsar rotational parameters in 
order to account for the effects of long-term timing noise. 

The pulsar timing observations were made primarily 
at frequencies around 1.4 GHz, employing a dual-channel 
cryogenic receiver system sensitive to two orthogonal cir- 
cular polarizations. The signals from each polarization 
were mixed to an intermediate frequency, fed through a 
multichannel filtcrbank, and digitized. The data were 
dedispersed in hardware and folded on-line according to 
the dispersion measure and topocentric period of the pul- 
sar. The folded pulse profiles were stored for subsequent 
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TABLE 1 
List of Calibrators 



Pulsar 


Calibrator 




R.A. 
(J2000) 




Dec. 
(J2000) 


Reference 


B0031- 


-07 


Tnn9zL— 041 9 


00: 




—04: 


■19-01 ^4874 


VCS 






Tn r^oa 
J-il Ucdlli 


00: 


■'^'^■97 777^^fi 


— 07: 


■19-'^^^ Q^^n9 
,±z.oo.you^ 




B0136- 


f57 


701 47+5840 


01: 


■47-4fi 541 1 49 


-1-58: 


■40-44 971 98 


VCS 






±11 Utidlli 


01: 


■'^7-'^n ZL7R9Q 


+58: 


■14-11 900^^ 




B0450- 


-18 


Tn4"S0— 1 837 


04: 




— 18: 


■37-00 40798 


VCS 






±11 Ucdlli 


04: 


■fi9-1 1 Q7QnS 
,u^.±±.y I yuo 


— 18: 


■1fi-41 Q^S9 




B0450- 


f55 


J U'-tOO^OOUO 


04: 




+55: 


OS-49 O'^SI 


VCS 


J0538+2817 


70^47-^2721 


05: 


■47-S4 148Q2S 


-I-27' 


■21 -^8 84282 


ICRF 






±11 Ucdlli 


05: 


.oo. ±u.O!joy 


+28: 


■94-9Q 4778 




B0818- 


-13 


TnS20— 1 258 


08: 


■20-'^7 447fi1fi 


— 12: 


■^8-'^Q lfiQ14 

.'jo.ijy.±'jy±rt 


ICRF 






±11 Utidlli 


08: 


■90-1 1 4'^^71 


— 13: 


■4Q-1 9 S'^'^8 




B1508- 


455 


J ±o ±Un^o I 


15: 


. ±\j .yj ^ . y Zi^ouo 




■09-4'^ '^7'^87 


ICRF 






Tri— rif>n tti 

±11 LJ*_.dlll 


15: 


■1 1 -48 0fi07fi 


+55: 


■41 ■'^fi 2876 




B1541- 


f09 


T1 f^c;f;4-i 1 1 1 

O l-OOO^ ± ± ± ± 


15: 




4-1 1 ■ 


1 1 -94 ■^fi'^71 


VCS 






Tri— rir>a in 

±11 LJt.dl±l 


15: 


■4'^-0'i 7021 4 


+09: 


■1 9-09 7845 




J17134 


-0747 


T1 71 q+n^l 7 


17: 


■1 Q-'^9 90fi91 f\ 


+08: 




VCS 






In-beam 


17: 


:13:27.09521 


+07: 


:30:18.9931 




B1933- 


H6 


J1924+1540 


19: 


:24:39.455870 


+15: 


;40:43.94172 


VCS 






In-beam 


19: 


:35:55.52636 


+16: 


:17:00.6849 




B2045- 


-16 


J2047-1639 


20: 


:47:19.667027 


-16: 


:39:05.84254 


VCS 






In-beam 


20: 


: 48:49.57191 


-16: 


:00:13.9195 




B2053+36 


J2052+3635 


20: 


:52:52.054980 


+36: 


:35:35.30036 


VCS 






In-beam 


20: 


:55:12.22739 


+36: 


:34:22.6929 




B2154- 


f40 


J2202-f4216 


22: 


:02:43.291372 


+42: 


;16:39.97992 


ICRF 






In-beam 


21: 


:57:54.02421 


+40: 


:04:54.6025 




B2310- 


1-42 


J2322+4445 


23: 


:22:20.358088 


+44: 


:45:42.35349 


VCS 






In-beam 


23: 


:12:13.49560 


+42: 


:39:32.3585 





Note. — Primary and in-beam calibrators observed for as trometry of pul- 
sars. The prim ary calibrators are either part of the ICRF I IMa et al.|[T99g: 
IFev et al.ll2004l) or tied to the ICRF by the VLB A Calibrator Surve y programs 
(VCSMBeaslev et al.ll2002l: IFomalont et al.ll200a :IPetrov et al.1120051 '1 . When in- 
beam sources arc listed, the pulsar positions are tied to the positions of the 
primary calibrator through the in-bcam sources. All in-beam calibrators were 
discovered as part of this project, as described in i|2] ICRF sources used here 
have positional uncertainties of ~ 200 ^tas. The sources from the VCS catalog 
have positions that are certain to 5 mas or better. 

etry. 



3.3. Correlation 

In order to mitigate time-average and bandwidth 
smearing (as for the VLA, ^2.1^ . the correlation was per- 
formed with 32 spectral channels of 250 kHz width and 
then averaged for 2 s. Both time-average and bandwidth 
smearing should be negligible over a 30' (FWHM) field of 
view; however, the ionosphere can induce unpredictable 
delays and rates that could otherwise reduce the correla- 
tion coefficients through decorrelation for greater degrees 
of averaging. Only the parallel senses of polarization 
were correlated. 

Each of the in-beam calibrators and pulsars were corre- 
lated in separate passes, with the phase reference center 
set in turn to the in-beam calibrator or pulsar positions. 
The pulsar correlation pass was gated as described above. 
Correlator parameters, such as clock offsets, station loca- 
tions, and Earth orientation parameters, were kept iden- 
tical to enable cross calibration between the correlator 
passes. 



4. VLBA DATA REDUCTION 

The multi-epoch observations of each pulsar typically 
produced 18 TB of raw voltage-sampled data that gener- 
ated 5 GB of correlated data, from which 5 astrometric 



analysis. 

Complete details regarding the procedure used to de- 
termine the pulsar rotational phases can be found in 
iHobbs et all (|2004f) : here we summarize only the salient 
aspects. Each timing observation consisted of a series of 
1-3 minute integrations, with a total duration of up to 
~45 minutes. The total observing time per pulsar was 
determined by the amount of time required to obtain 
sufficient sensitivity, given the pulsar flux density. In 
post-processing, any of the individual integrations that 
appeared to be dominated by RFI were removed, the po- 
larizations combined, and the data were averaged to pro- 
duce a single total-intensity profile for the observation. 
Pulse times of arrival (TOAs) were subsequently deter- 
mined by finding the peak in a time domain convolution 
of the averaged profile with a pulse profile template corre- 
sponding to the observing frequency. These TOAs were 
corrected to the solar system barycenter using the Jet 
Propulsion Lab oratory DE200 solar system ephemeris 
(|Standish|[l983) . analyzed using the standard TEMPO 
software package^^, and transmitted to the VLBA corre- 
lator for use in the correlation. We note that our refined 
astrometry does not affect the pulse phase prediction 
from the timing solutions at a relevant level, and while 
gating improves the S/N ratio of the pulsar in VLBA 
observations, it does not otherwise influence the astrom- 

|http: //www. atnf . cslro . au/research/pulsar/tempo/ | 
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parameters {ao, Sq, ^a, fJ-s, tt)^"^ are finally derived; the 
factor of ~ 10^^ truly warrants the use of the term "data 
reduction" for VLBI pulsar astrometry. Given the data 
volume, a semi-automated pipeline was built within the 
NRAO Astronomical Image Processing System (AIPS^^) 
for calibration and imaging. The pipeline allowed for 
human interaction during data editing and specific imag- 
ing steps while minimizing errors during repetitive book- 
keeping and calibration steps. We summarize the vari- 
ous calibration steps below, with emphasis on the specific 
corrections introduced. 



4.1. A priori Calibration 

Corrections for various instrumental and system pa- 
rameters are determined from geometric factors and an- 
cillary data. Since they do not require observations of 
an astronomical source, we distinguish them as a priori 
calibration, following usual practice. These corrections 
are listed below, in the order in which they were applied. 

Earth Orientation Parameters: In order to combine 
the signals from the various antennas, the correlator 
must have an accurate geometrical model of the array 

omnevlfl999t ). The model includes the orientation of 
the Earth in space, as described by a set of Earth orien- 
tation parameters (EOPs). The most accurate estimates 
of these EOPs are typically available only after a time lag 
of several weeks, well after data correlation. Moreover, 
during the course of this project, it was discovered that 
the VLBA correlator had inadvertently been provided 
with predicted rather than measured EOPs (Walker et 
al. 2005, VLBA Test Memoes 69). The resultant dif- 
ferential astrometric errors are potentially as large as 
~ 2 mas, dwarfing the astrometric error budget. The 
correction for this was a delay offset computed between 
the delay model using the (incorrect) correlation EOPs 
and the actual EOPs, and this correction was applied to 
each data set. 

Antenna Positions: The locations of the VLBA anten- 
nas have an impact on absolute astrometry at the level 
of ~ 50 /ias per mm of position error on the longest base- 
lines; the astrometric consequences are worse on shorter 
ones. The errors in differential astrometry are smaller 
than this by roughly the target-calibrator separation (in 
radians) or ~ 1 /zas per mm of position error per degree of 
separation on the sky. On 2004 May 12, a new model for 
antenna positions and velocities (GSFC solution 2004b) 
was adopted at the VLBA correlator, introducing a typi- 
cal discontinuity of ~ 20 mm in antenna positions. Each 
data set correlated before the change was corrected by 
shifting antenna positions to be consistent with the 2004b 
solution. 

Quantization Corrections: Because the VLBA uses 
two- or four-level quantization when sampling, the visi- 
bilities produced by the VLBA correlator differ system- 
atically from those that would be produced by an ideal 
analog system. Quantization statistics were calculated 
for the auto-correlations and used to correct the ampli- 
tudes of the cross-correlations. 

^•^ In this work, both fia and fig always refer to angular, rather 
than coordinate, rates. Specifically, fia = ^ cos S. 
http: //www. aips .nrao ■ edu 
^® See |http : //www ■ vlba . nrao ■ edu/memos/test / 1 



Gain and System Temperature: The measured system 
temperatures and antenna gains (typically ^0.1 K Jy~^ 
for a VLBA antenna at 1.5 GHz) were used to scale the 
correlation amplitudes, in the standard manner. 

Parallactic Angle Correction: The VLBA antennas 
have altitude-azimuth mounts, so the antenna feeds ro- 
tate as the antennas track sources across the sky. Cor- 
rections were applied for the geometric phase shift in- 
troduced between the left and right circularly polarized 
signals due to feed rotation. 

Coarse Ionospheric Correction: Precision astrometry 
requires correction for the differential phase and delay 
effects introduced by signal propagation through the 
Earth's ionosphere. The effects due to different iono- 
spheric conditions above each antenna fall into two cat- 
egories: image distortion and image shift. The former 
causes an overall decrease in the signal-to-noise of the as- 
trometry by broadening the image. The latter introduces 
an unknown, unwanted position offset, directly impact- 
ing the astrometric results. 

The ionospheric total electron content (TEC) can be 
estimated by measuring the propagation delay between 
Global Positioning System (GPS) signals at two different 
frequencies along the lines of sight to the constellation of 
GPS satellites. Global maps of the ionospheric TEC are 
routinely produced based on such measurements from 
around the world, and archived by the NASA Crustal 
Dynamics Data Information System (CDDIS^^). These 
maps are relatively coarse, typically quantifying the iono- 
spheric TEC on a grid with pixels spaced 5° in longi- 
tude and 2? 5 in latitude, and updated every 2 hr. While 
these coarse models are far from sufficient for precision 
astrometry, Walker & Chatterjee (1999, VLBA Scientific 
Memo^'^ 23) use them to find a reduction in phase delay 
scatter by a factor of 2 to 5 at 2.2 GHz, and using them 
to correct for ionospheric dispersion delay did improve 
phase referencing. 

4.2. Radio Frequency Interference Identification and 
Excision 

Even though the spectral windows were chosen to min- 
imize the overlap with known sources of interference, 
some residual RFI persisted in the data. While the usual 
interferometric practice is to edit data on a per-baseline 
basis, we concluded that RFI could be identified on an 
antenna basis, i.e., identified with being present only on 
baselines common to a specific antenna. Therefore, we 
could identify and flag RFI based on the autocorrelation 
data at each antenna, and baseline-based RFI excision 
was only rarely needed. 

Data inspection and editing required extensive human 
intervention. However, since the same data were corre- 
lated multiple times, we could identify and flag RFI on 
the output of the flrst correlation pass and apply the flags 
to visibilities generated on all further passes. 

4.3. Astronomical Calibration 

Phase, delay, and bandpass calibrations for VLBI are 
typically derived from observations of calibrator sources. 
For a given pulsar the same bright VLBA calibrator, as 
listed in Table [H was used for all three purposes, and we 



' |http : //cddisa . gsf c ■ nasa ■ gov7 
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followed an iterative process of calibration refinement. 
The calibration steps listed below were all performed us- 
ing MPS. 

Fringe fitting: We solved for calibrator visibility phase 
delays and time derivatives (rates) by fitting for a phase 
slope across each frequency band, assuming a point 
source model. The delays were smoothed by averaging 
over 15 minutes and interpolated for observations of the 
targets. 

Bandpass calibration: After fringe-fitting, calibrator 
data from all epochs on a particular target were com- 
bined, and the inner 26 of the 32 channels in each fre- 
quency band were then imaged. Several iterations of 
imaging and self-calibration were performed to create 
a final calibrator image, averaged over all observation 
epochs. The resulting 'global' calibrator model was used 
for bandpass calibration at each epoch, determining a 
phase and amplitude gain factor for each of the 32 chan- 
nels of each frequency band, for each antenna. These 
corrections help to minimize baseline-dependent errors 
in further calibration steps. 

Frequency division: The calibrator structure did not 
change significantly over the duration of the project, en- 
abling us to build time-averaged models. However, in 
some cases the calibrator images showed significant fre- 
quency dependent structure. As such, after the initial 
fringe-fitting and bandpass calibration, we used separate 
calibrator models for each frequency band. 

Phase and amplitude calibration: For each frequency 
band, calibrator data from all epochs were combined, and 
an imaging and self-calibration loop was used to derive 
an epoch-averaged calibrator model. The global calibra- 
tor model produced for bandpass calibration was used as 
the initial model in each case. Phase and amplitude cali- 
brations were then derived independently for each epoch 
and frequency band, using calibrator models that were 
averaged over all epochs but not over frequencies. The 
phase and amplitude calibration was then applied to the 
target source (i.e., either the pulsar or the in-beam cal- 
ibrator) and written to a single-source, calibrated UV 
database. In order to ensure that the data for the pul- 
sar and in-beam calibrator(s) were identically calibrated, 
we derived the fringe-fitting, bandpass calibration, and 
phase/amplitude calibration for one correlation pass and 
copied these tables to all subsequent correlation passes. 

At the end of this calibration sequence there were typi- 
cally 32 calibrated single-source UV databases per pulsar 
and per in-beam calibrator, one for each of 4 frequency 
bands at 8 epochs. 

4.4. In-beam calibration 

The angular separation between an in-beam calibrator 
and the target pulsar (typically ^^15-20') is much smaller 
than their angular separation from the phase calibrator 
(typically ~2-3°). As such, most of the residual calibra- 
tion errors will be common to the pulsar and the in-beam 
source, and we can exploit in-beam calibrators to incre- 
mentally improve the calibration of the pulsar. 

As for the phase calibrators, we created indepen- 
dent models for in-beam calibrators at each frequency 
band by combining and imaging the calibrated visibil- 
ity data from all epochs at that band. A few passes of 
self-calibration (usually phase-only, but including ampli- 
tude self-calibration for in-beam sources brighter than 



50 mjy) were used to improve the models. Since in- 
beam calibrators are usually fairly weak (~ 10 mJy) and 
show some source structure, we required longer time av- 
erages and lower S/N cutoffs for the calibration solutions, 
up to a limit of 5 minutes of averaging and S/N > 1.8. 
This lower S/N threshold was determimed to maximize 
the S /N in the phase-referenced image of the pulsar af- 
ter the phase transfer from the in-beam calibration so- 
lutions; calibrators with different characteristics may be 
optimized with slightly different cutoff values. In some 
cases, the in-beam source had enough resolved structure 
that the longest baselines of the VLBA (to Saint Croix 
and Mauna Kea) had to be discarded. Even in such cases, 
using in-beam calibration provided superior astrometric 
results compared to retaining the longest baselines while 
only referencing observations to the more distant phase 
calibrator. 

With an epoch-averaged model for the in-beam cali- 
brator at each frequency band, we used a final pass of 
phase-only self-calibration to obtain phase corrections. 
These corrections were applied to the calibrated pulsar 
data and used for subsequent imaging. 

4.5. Imaging and position fitting 

Since the target pulsars are intrinsically compact, the 
final imaging step is quite simple. The VLBA beam 
is North-South elongated, with typical dimensions of 
5 X 14 mas for our chosen frequencies, observation dura- 
tions, and selection of source declinations. We used pixel 
sizes of 1 X 1 mas and approximately natural weighting 
(an AIPS 'robustness' parameter of 2, where 5 corre- 
sponds to pure natural weight ing and —5 corresp onds to 
uniform visibility weights; see lBriggs et al.|[l999( ). Many 
of the target pulsars have observed diffractive scintilla- 
tion timescales at 1.5 GHz that are comparable to the 
2 hr observation span, and so their sidelobes are often not 
well subtracted. For this reason support constraints (i.e. 
'boxes') were used to restrict the region to be cleaned, 
and we stopped the cleaning process interactively once 
the residuals within the support box were at the level of 
the image noise far from the pulsar. Finally, we deter- 
mined astrometric positions for our targets at each epoch 
and frequency band by fitting elliptical Gaussians in the 
image domain. 

5. ANALYSIS AND RESULTS 

At the end of the calibration and imaging process, we 
are left with a set of positions {at,!/, 5ty, t = 1 . . v — 
1 ... 4} at each epoch t and frequency band v. We fit 
these positions with a five-parameter astrometric model 
encompassing position (aoi'^o), proper motion (/iQ,,/!^), 
and parallax (tt) using a standard weighted linear least 
squares technique. The key problem is the determination 
of the weights, and hence the fit uncertainties, in the 
presence of unknown systematic errors at each epoch. 

At each epoch, the target position errors have a (ran- 
dom) measurement error component, corresponding to 
the beam FWHM/(2 x S/N), where S/N is the signal- 
to-noise ratio of the final pulsar image. However, using 
only these random errors produces fits with reduced 
values > 1, implying the presence of unmodeled system- 
atic errors. In VLBI astromctry at 1.5 GHz, the biggest 
systematic error component is usually due to the resid- 
ual unmodeled ionosphere between the in-beam calibra- 
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tor and the target pulsar, with other possible contribu- 
tions from errors in the Earth orientation parameters 
(c.f. ii4.ip . calibrator position errors, or structural evo- 
lution of the calibrators. The calibrator position errors 
could in principle impart a significant ('^ 50 /las) time- 
dependent error on the phase-referenced position mea- 
surements. (We note that uniformly observing a given 
target over the same hour angle range causes calibrator 
source position errors to produce a consistent astrometric 
offset at each epoch, affecting only the absolute position 
determination.) A common approach to treating unmod- 
eled errors is to inflate all the random errors by a trial 
factor until a reduced ^ value f« 1 is obtained. However, 
the residual ionospheric errors, for example, depend on 
the proximity of the Sun (c.f. i j3.ip . and since the sep- 
aration varies with season, scaling the random errors at 
each epoch by the same factor would not be optimal. In 
cases where the angular separation between the sun and 
the target is more than 30°, the uncalibrated ionosphere 
typically yields astrometric errors of ^ 200 /^as for a 1° 
target-calibrator separation at an observing frequency of 

1.5 GHz. 

In our past programs (e.g., IChatterjee et al.l l2001l . 
l2004f) we have investigated several possible measures of 
the epoch-dependent systematic error, including the an- 
gular size of the pulsar image and its deviations from the 
beam shape, or the scatter between the positions deter- 
mined at each frequency band at a given epoch. While 
these measures are useful, they are not perfect, and re- 
gression analysis between the estimated systematic error 
and parameters such as the Sun angle or the elevation 
of the source above the horizon do not reveal a simple 
correlation. 

Instead we choose to estimate the fit uncertainties 
from the position measu rements th emselves using the 
bootstrap technique (e.g.. lEfron fc T ibshirani 1991il). We 
treat the iV^ data points \a.t,v,^t,v] at each epoch and 
frequency band as independent estimates of the (time 
dependent) position, and select samples from the 
set with replacement, discarding degenerate samples with 
fewer than three independent positions. We then fit 
the sample for astrometric parameters using standard 
weighted least squares, using only the random measure- 
ment uncertainties. The process of sampling and fitting 
is repeated a large number of times. Each resulting fit 
generally has a reduced > 1, but only the values of 
the fit parameters are of interest. Histograms of the re- 
sulting fit parameters are shown for a typical case with 
50,000 sampling iterations (PSR B2310-I-42, Figure [l]). 
It is apparent that the parameters follow a smooth dis- 
tribution with a well-defined most probable value (the 
distribution mode). The most probable value for each 
astrometric parameter and their most compact 68% con- 
fidence intervals as estimated from the bootstrap method 
are listed in Table [2] for the fourteen pulsars in our final 
sample. The most probable values and confidence inter- 
vals for the distances and transverse velocities of these 
pulsars were also derived directly via bootstrap, and are 
listed in Table [31 

In all cases, the values of the parameters estimated by 
the bootstrap method were consistent with the param- 
eters obtained by simple least squares fitting (using all 
epochs without replacement) to well within the bootstrap 
68% confidence interval. However, the uncertainties in 




24 24.1 24.2 24.3 5.6 5.8 6 6.2 0.6 0.9 1 1.1 
Me [mas/yr] [mas/yr] n[mas] 



Fig. 2. — Uncertainties in astrometric parameters as estimated 
by the bootstrap method for PSR B2310-(-42. The top row shows 
histograms of the estimates for the fit parameters from 50,000 boot- 
strap iterations, equivalent to an un-normaUzed probability distri- 
bution function (PDF), and the vertical line indicates the most 
common value (the mode). The bottom row shows the cumulative 
distribution function (CDF) derived from the histogram above, and 
the vertical line indicates the median parameter value. In each 
panel, the horizontal bar indicates the most compact 68.26% confi- 
dence interval (i.e., the smallest interval that includes 34,130 esti- 
mates for each given parameter). In this and every other case, the 
median and mode are consistent to well within the 68% confidence 
interval. The proper motion in right ascension, /la , is the angular 
rate, ^ cos 5. 

' at 
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Fig. 3. — Uncertainties in astrometric parameters as estimated 
by the bootstrap method for PSR J1713-I-0747. As in Figure |2] 
the top row shows the histogram of the estimates (equivalent to 
an un-normalized PDF) and the bottom row shows the cumula- 
tive distribution function (CDF) for each astrometric parameter, 
while vertical lines indicate the mode and median estimates and 
horizontal bars show the most compact 68% confidence interval for 
each parameter. The fit for fia is the only case among all pulsars 
discussed in the present work where a bimodal distribution was 
obtained for an astrometric parameter, and thus represents our 
worst-case scenario. However, our astrometric results for this pul- 
sar can be compared to independent estimates from pulse timing: 
see i|6.3l and Table |4] for a discussion. 

the estimates of the astrometric parameters differed be- 
tween the two methods. The power of the technique is 
illustrated by our worst-case scenario (PSR J1713+0747, 
Figure 121), which was the only case where a bimodal dis- 
tribution was obtained for any of the parameters. The 
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bimodality in /Iq, reflects the presence or absence of one 
well-measured position, and causes the uncertainties to 
deviate significantly from symmetry. However, there is 
no ambiguity in the most probable parameter values or 
in the bootstrap confidence intervals. Coincidentally, 
PSR J1713-I-0747 has well determined astrometry from 
pulse timing as well, allowing a stringent check on our 
results, as we discuss further below (i i|6.3|) . 

Compared to simple least squares fitting, the boot- 
strap method is expected to return larger fit uncertain- 
ties, since the time sampling for each trial in the boot- 
strap process is generally worse than simply using all 
epochs of data in the fit process. To compare the perfor- 
mance of the bootstrap method to simple least squares 
fitting, we simulated observations of PSR J1713-f 0747 on 
the actual observation dates, and recovered astrometric 
parameters both by a least squares fit to all the data, 
and by bootstrap. As expected, for Gaussian random 
uncertainties in the input position at each epoch, the 
two methods returned astrometric parameters that were 
essentially identical, differing by much less than the fit 
uncertainties, and the bootstrap method returned some- 
what larger confidence intervals compared to the simple 
least squares fit (between and 25%). Next, we added 
unmodeled systematic errors to the simulated positions, 
such that the mean "observed" position at a given epoch 
differed from the "true" position and the observation un- 
certainties understated the effective simulated errors. As 
systematic errors were added to increasing numbers of 
epochs, the results of the two techniques began to di- 
verge. While the recovered best-fit astrometric parame- 
ters did not differ significantly between the two methods, 
the bootstrap fit uncertainty intervals generally encom- 
passed the "true" parameter values (within 1- or 2-a) 
while the simple least squares method returned uncer- 
tainties that were too optimistic and failed to include 
the input values within several multiples of a with in- 
creasing systematic errors. These simulations therefore 
justify our decision to use the bootstrap method when 
we are not able to reliably quantify the systematic errors 
in our observations; while our astrometry may be some- 
what less precise as a result, it is more likely to be more 
accurate. 

For illustration purposes, we show the parallax and 
proper motion signature of PSR B1933+16 in Figure [H 
and the parallax signature of PSRs B0818— 13 and 
B2310-f42 in FigureO In these figures, the position error 
bars include estimates of the systematic error from the 
scatter between positions measured at each epoch at dif- 
ferent frequency bands. However, we note that these es- 
timates are used only for illustrative purposes, and the fit 
uncertainties are obtained from bootstrap, as described 
above. 

Finally, we note that signal-to-noise considerations 
required the averaging of all frequency bands for 
PSR B0450— 18. As such, the bootstrap was performed 
with only 8 position measurement pairs {at , St at each 
epoch) while fitting for 5 astrometric parameters, re- 
sulting in a high proportion of degenerate samples. For 
this one object, we repeated the bootstrap for 10^ iter- 
ations, but obtained no difference in the results or their 
confidence intervals compared to 10^ iterations. How- 
ever, the uncertainties for the astrometric parameters for 
PSR B0450— 18 are much larger than for the rest of our 
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Fig. 4.— The proper motion and parallax of PSR B1933+16. 
The best fit parallax and proper motion (fia = 1.1 mas yr~^, 
fig = —16.1 mas yr~^, tt = 0.2 mas) is overplotted. The earliest 
position (top left) is at MJD 52485, while the final position (bottom 
right) is at MJD 53116. ARA is the angular displacement Aa cos 5. 
Note the different scales on the two axes. 

sample. 

6. DISCUSSION 

We have presented parallaxes and proper motions for 
14 pulsars (Table [2]) and derived estimates for their dis- 
tances and transverse velocities (Table [3]). These results 
represent a significant increase in the number of pulsars 
with trigonometric parallaxes, and include distances to 
two pulsars over 5 kpc away that are determined to bet- 
ter than ~15% : most probable D = 7.2~^\'^ kpc for 

PSR B154H-09, and 5.5tJ:g kpc for PSR B2053+36. 
Early results from this program were published for two 
objec ts in the sample, PSR B 15 08-1-55 (IChatteriee et all 
[2005h and PSR J0538-f2817 ( Ng et al.l 120071 ). but with 
a simpler, more ad hoc treatment of the systematic er- 
rors. The bootstrap approach followed here does not as- 
sume that errors at each epoch are normally distributed, 
a condition that is now appreciated to be rarely true in 
phase referenced astrometry at frequencies below 5 GHz, 
where the turbulent ionosphere dominates the phase- 
referencing errors. For situations with over-determined 
measurements, such as those presented here, the boot- 
strap method effectively results in discrepant data re- 
ceiving a lower weight, and so the current results are 
more robust. 

For PSR B1508-I-55, the astrometric parameters de- 
termined here are consistent with our previous work, 
and the potential birth site in the Cygnus OB associ- 
ations is unchanged. Nor is the high estimated veloc- 
ity {V± « 1000 km s-i; Vsd > HOO km s^^) signifi- 
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TABLE 2 

ASTROMETRY RESULTS 



Pulsar 


RA 


Dec 


Ma M(5 


7r 




(J2000) 


(J2000) 


(mas yr"-"-) (mas yr"-"-) 


(mas) 



B0031- 


-07 


uu 




UO 


O ( uot^ i ) 


07 


91 


oo 




1 n 


,7+0.08 
■^'-0.08 


1 1 
— ±± 


,0+0.11 
^'-'-0.16 


u 


B0136- 


f57 


Ul 


oy 


1 Q 


/ ^un^ iz j 




1 A 


o± 


o±y i / ) 


1 Q 

— ly 


1 1 +0.06 
-^-^-0.07 


1 f\ 
— ±o 


fin+o-os 


n 
u 


B0450- 


-18 


04 


52 


34 


1057(1) 


— 17 


59 


23 


371(2) 


8 


"'^ — 2.24 


10 


"■^-1.29 





B0450- 


f55 


04 


54 


07 


7506(1) 


-f55 


43 


41 


437(2) 


53 


0^+0.06 
■^^-0.05 


-17 


trfj+0.12 
'^'-'-0.14 





J05384 


-2817 


05 


38 


25 


0572(1) 


-1-28 


17 


09 


161(2) 


-23 


£•7+0. 10 
O'-O.IO 


52 


O7+0.09 
°'-0.10 





B0818- 


-13 


08 


20 


26 


3817(1) 


-13 


50 


55 


859(2) 


21 


64+0-09 
"^-0.09 


-39 


44+0-04 
^*-0.05 





B1508- 


f55 


15 


09 


25 


6298(1) 


-1-55 


31 


32 


394(2) 


-73 


64+0-05 
"^-0.04 


-62 


fic:+0.09 
"3-0.08 





B1541- 


f09 


15 


43 


38 


8250(1) 


-f09 


29 


16 


339(2) 


-7 


CI +0.06 
"-■--0.05 


-2 


07+0. 06 
°'-0.07 





J17134 


-0747 


17 


13 


49 


5306(1) 


-1-07 


47 


37 


519(2) 


4 


7t;+0.17 
'3-0.07 


-3 


f.7+0.16 
-0.15 





B1933+16 


19 


35 


47 


8259(1) 


-fl6 


16 


39 


986(2) 


1 


,0+0.12 
^■^-0.13 


-16 


(,0+0-15 
^^-0.13 





B2045- 


-16 


20 


48 


35 


6013(2) 


-16 


16 


44 


534(3) 


113 


1O+0.02 
^"-0.02 


-4 


(,(,+0.28 


1 


B2053- 


f36 


20 


55 


31 


3521(1) 


-f36 


30 


21 


469(2) 


1 


04+0-04 
^*-0.04 


-2 


46+0-13 
*"-0.13 





B2154- 


f40 


21 


57 


01 


8495(1) 


4-40 


17 


45 


986(2) 


16 


1 O+0.09 
^■^-0.10 


4 


-,9+0.12 
^^-0.12 





B2310- 


1-42 


23 


13 


08 


6209(1) 


-1-42 


53 


13 


043(2) 


24 


1 t;+0.10 
-'-'J-0.07 


5 


Qt:+0.13 
='3_o.l2 






1+0.08 
^-0.07 
,+0.04 
'-0.04 
- + 1.40 
^-0.59 
1+0.04 
*-0.05 
,+0.12 
^-0.09 
1+0.03 
"^-0.04 
7+0.03 
'-0.03 
,+0.02 
^-0.02 
-+0.06 
^-0.05 
-1+0.08 
^-0.12 
-+0.03 
^-0.02 
7+0.03 
'-0.03 
S+0.06 
'-0.06 
,+0.06 
^-0.07 



Note. — Pulsar coordinates have been adjusted for improvements in primary calibrator positions. 
Uncertainties in these coordinates (in parentheses) reflect uncertainties in the positions of primary 
calibrator; a minimum uncertainty of 0.0001 second of time in right ascension and 2 mas in declination 
are assigned as calibrator source structure evolution with frequency could cause systematic offsets from 
the reported 8.4 GHz positions at about this level. In all cases this calibrator position uncertainty is 
thought to be greater than the relative error introduced in phase-referencing. Quoted proper motion 
and parallax uncertainties are most compact 68% confidence intervals. All astrometry is performed 
and reported in Equinox J2000. Reported positions are for epoch 2002.0. 



TABLE 3 

Parameters derived from astrometr.y 



Pulsar 


/ 


b 


DM 


Ddm 






Vx 




(°) 


n 


(pc cm~^) 


(kpc) 




(kpc) 


(km s-i) 


B0031-07 


110.42 


-69.82 


11.38 


0.4 


1 


^^"-0.09 


77l« 


B0136+57 


129.22 


-4.04 


73.78 


2.8 


2 


fic:+0.35 


319lg 


B0450-18 


217.08 


-34.09 


39.90 


2.4 





7C+0.36 
'"-0.60 


54lg 


B0450+55 


152.62 


7.55 


14.50 


0.7 


1 


iq+0.07 
^^-0.06 


31711° 


J0538+2817 


179.72 


-1.69 


39.57 


1.2 


1 


on+0.22 
■^"-0.16 


357+5° 


B0818-13 


235.89 


12.59 


40.94 


2.0 


1 


06+017 
y"-0.12 


418li 


B1508+55 


91.33 


52.29 


19.61 


1.0 


2 


ir,+0.13 
J^^-0.14 


9631^1 


B1541+09 


17.81 


45.78 


35.24 


>35 




7.211-3 


077+52 

^"-42 


J1713+0747 


28.75 


25.22 


15.99 


0.9 


1 


n^+o.de 

""-0.07 


30+2 


B1933+16 


52.44 


-2.09 


158.52 


5.6 




"■^-2.7 


004+108 
"^^-208 


B2045-16 


30.51 


-33.08 


11.46 


0.6 





Qc:+0.02 
''"-0.02 


512+1^ 


B2053+36 


79.13 


-5.59 


97.31 


4.6 




5 5+1-2 
"■"-0.8 


6911^ 


B2154+40 


90.49 


-11.34 


70.86 


3.7 




4+0.9 
"■^-0.7 


2641^3 


B2310+42 


104.41 


-16.42 


17.28 


1.2 


1 


(,(,+0.08 
^'"-0.07 


1251^° 



Note. — Pulsar distances DjjM were estimated fro m their disp ersion mea- 
sure (DM) using the NE2001 electron density model HCordes fc L azio 20021). 
Confidence intervals for the parallax distances Dtt and transverse velocities V± 
were estimated directly from the bootstrap. 



cantly altered. For PSR J0538+2817, the proper mo- 
tion in declination received an adjustment of about 2 a 
when employing the bootstrap analysis, but the magni- 
tude of the change is small (from 52.59 ± 0.13 mas yr~^ 
to 52.87™ mas yr and does not affect the conclu- 
sions of iNg et all (|2007t l. The re sults also remain con- 
sistent with the timing results of iKramer et al.l (|2003l ). 
but are much more precise. Below we discuss some of 
the implications of our new measurements, particularly 
for birth sites, kinematic ages, Galactic electron density 



models, and reference frame ties. 

6.1. Pulsar Birth Sites and Kinematic Ages 

Accurate astrometric measurements allow pulsars to 
be trace d back through the Galaxy to potential birth 
sites (e.g. iHoQgerwerf et al.ll200Cll[200lHVlemmings et alj 
l2004f ). If a birth site can be identified, the three- 
dimensional birth velocity can be estimated, a kinematic 
a-ge Tkin can be inferred, and a combination of the birth 
spin period and braking index can be constrained. A 
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Fig. 5. — The parallcix signature of two pulsars in Right Ascension and Declination, after subtracting the respective best-fit proper motions 
from the astrometric positions. ARA is the angular displacement Aa cos (5. Left; The parallax of PSR B0818— 13, with curves overplotted 
corresponding to the best fit parallax tt = 0.51 mas. Right: The parallax of PSR B2310+42, with curves overplotted corresponding to the 
best fit parallax tt = 0.93 mas. In each panel, the inner error bars indicate the random position uncertainties, while the outer error bars 
indicate the net estimated uncertainties (random and systematic, added in quadrature) for illustration purposes only. 



detailed analysis of the Galactic trajectories of all pul- 
sars with precise astrometry is beyond the scope of this 
work and will be presented in future (Vlemmings et al., 
in prep.), but a preliminary assessment of the current 
sample of pulsars is presented here. 

The trajectories of the pulsars in our sample (except 
the old, recycled pulsar J1713-I-0747) are shown in Fig- 
ure ini for three different values of the (unknown) radial 
velocity {V, = -200,0 and +200 km s'^). The pulsar 
trajectories were calcu lated using the t hree comp onent 
Stackel potential from lFamaev &: Deionghd (I2003D with 
the sc aling parameters as described in iVlemmings et aLl 
()2004[ ). Figure [S] also shows the po sitions and ex- 
tent o f the nearby OB associations from Ide Zeeuw et al.l 
(|2_999l) and a number of open clusters compiled from the 
WEBDA catalog^* with a distance < 5 kpc and a well- 
defined age < 10 Myr. We only plot those open clusters 
for which a pulsar trajectory crosses within 20 pc of the 
cluster center on the sky, without accounting for the dis- 
tance or the motion of the cluster itself. The motion of 
the cluster can be several hundred pc during the life- 
time of a pulsar, and needs to be taken into account in 
a more detailed analysis. Although a number of pulsar 
trajectories pass close to the OB associations or open 
clusters on the sky, the pulsar distance is typically in- 
compatible with an origin in the cluster or association. 
However, besides the earlier identification of the likely 
birth site of PSR B1508-F55 in one of the Cygnus OB 
associations (jChatteriee et al.l 120051) , the coarse analy- 
sis performed here reveals the possible birth location of 
PSR B2045-16 in the open cluster NGC 6604, as dis- 
cussed below. 

The progenitors of NS are massive O- and B-stars, 
which have a low Galactic scale height ^63 pc, much 
smaller than the characteristic distance of radio pul- 
sars from the Galactic plane. This is consistent with 
pulsars having significantly larger peculiar velocities 

http:/ /www.univie.ac.at/webda/ 



than their proge nitors, as recognized for example by 
iGott et all (|l970[ ). As such, it is expected that most 
(young) pulsars shoul d be moving awa,y from the Galac- 
tic plane. Although ' Harrison et al.l |l99§) found that 
17% of the young pulsars in their sample seemed to 
be moving towards the plane, later analysis indicates 
that after accounti ng for the pulsar bh'th s cale height 
(zo = 160 ± 40 pc, lArzoumanian et al.|[200^ . all young 
pulsars do see m to be moving away fro m the Galac- 
tic plane (e.g.. |Cord cs & Chernoff 1998; iBrisken et al.l 
l2003at iHobbs et al.ll2005l) . Any apparent motion toward 
the Galactic plane can often be explained by the un- 
known radial velocity Vi, as the velocity perpendicular 
to the plane T4 is given by 

V; = Vb cos b + Vr sin b, (1) 

where Vb is the component of the velocity in the tangent 
plane along Galactic latitude b. 

Alternatively, pulsars older than ~ 20 Myr could po- 
tentially already be falling back towards the plane due 
to the effect of the Galactic potential on their trajectory. 
Figure [7] indicates against the location with respect 
to the Galactic plane for the pulsars in our sample. It is 
apparent that for a reasonable range of 14, all the pulsars 
are moving away from the Galactic plane, with the ex- 
ception of PSR J05384-2817, which is located well within 
the pulsar birth scale height a.t z — —41 pc. Figure[7]also 
shows that, as expected, the youngest pulsars are found 
closest to the plane. 

Assuming that all pulsars are indeed born in the Galac- 
tic plane between —160 < zq < 160 pc, we can estimate 
the age of the pulsar from the measured z = D sin b and 
Vz- The kinematic age Tkin is given by: 

mn = iDsmb-zo)/V,. (2) 

The estimated kinematic ages of the pulsars in our sam- 
ple are shown in Figure[51 for —200 < Vi- < -1-200 km s^^. 
As deceleration in the Galactic potential will reduce T4 
as pulsars get older, Tkin will be a lower limit on the 
true age. Spin-down ages Tsd ~ P/2P are also shown 
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Fig. 6. — The distribution of our pulsar sample in Galactic 
coordinates. The solid dots denote the present pulsar positions, 
and the lines represent their trajectory on the sky calculated for 
an age corresponding to the pulsar spin-down age. The open stars 
indicate the possible birth places of the pulsars. The three panels 
a,b and c represent the calculated orbits using Vr = +200,0 and 
—200 km s~^ respectively. From left to right (decreasing / of 
the present position) the pulsars are J0538+2817, B0450+55, 
B0136+57, B0031-07, B2310+42, B1508+55, B2154+40, 
B2053+36, B1933+16, B2045-16, B1541+09, B0818-13 and 
B0450— 18. The shade d areas indicate the nearby OB associations 
\Ae Zeeuw et al.1ll999I V From left to right and from top to bottom: 
Per OB2, a Per, Cep OB3, Cep OB6, Cep OB2, Lac OBI, 
Cyg 0B7, Cyg OB4, Set OB2, Upper Sco, Upper Cen Lupus, 
Lower Cen Crux, Tr 10, Vela OB2, Col 121 and Ori OBI. 
The open circles denote the open clusters from the WEBDA 
catalog with well-defined distances (< 5 kpc), positions and ages 
(< 50 Myr). Only those open clusters are shown for which a 
pulsar orbit approaches the cluster position within < 20 pc on the 
sky. 

for comparison in Figure [51 where P is the period of 
the pulsar, assumed to have been born spinning rapidly 
and slowing steadily with a braking index n = 3. Each of 
these assumptions (rapid spin at birth, steady slowdown, 
n = 3) is known to be problematic iii speci fic cases (e.g., 
iGaensler fc Frailll2000l : lKramer et "alll2003f) . However, in 
spite of the variety of assumptions built into both Tgd and 
Tkiii) the agreement between them for most of the pulsars 
in our sample is remarkable. 

The only pulsars for which there is an ap- 
parent discrepancy between the two age estimates 
are PSRs J1713-h0747, B0818-13, J0538-H2817 and 
B0450— 18. These pulsars appear to be falling towards 
the Galactic plane and thus have undetermined Tkin for 
K- = 0. The recycled pulsar J1713+0747 has hkely al- 
ready completed one or more orbits through the Galac- 
tic plane, and the simple analysis outlined above is 




[km/s] 

Fig. 7. — Pulsar distance above the Galactic plane (z) versus 
its velocity perpendicular to the plane (Vz)- The different symbols 
represent different ranges of pulsar spin-down age and the open 
star indicates the recycled pulsar J1713-I-0747. As 14 depends on 
the unknown radial velocity of the pulsar, the bars denote —200 < 
Vr < -1-200 km s~^. The pulsars in the shaded areas are moving 
away from the Galactic plane while those in the open areas are 
falling towards the plane. The horizon tal dashed lines indicate th e 
pulsar birth scale height determined bv lArzoumanian et al.l (120021 ). 
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Fig. 8. — Kinematic ages for our pulsar sample. The solid sym- 
bols denote the kinematic ages calculated using the pulsar distance 
above the Galactic plane and velocity 14, assuming Vr = km . 
The bars indicate the range -200 < Vr < -1-200 km s"!. The open 
circles are the spin-down ages. For those pulsars without solid 
symbols, the perpendicular velocity is pointing towards the plane 
when Vr = km s~^. As deceleration in the Galactic potential will 
have affected Vz for the pulsars older than ~ 10 Myr, ritin will be 
a lower limit. 

obviously inapplicable in that case. Allowing a ra- 
dial velocity slightly larger than the arbitrary cutoff of 
200 km s-i resolves the discrepancy for PSR B0818-13. 
PSR J0538-h2817 is still within the Galactic pulsar birth 
scale height, as indicated above. It is associated with the 
known supernova remnant S147, and the pulsar age of ^ 
30 kyr derived from the association (.Kramer et al.|[2003l ) 
is also quite inconsistent with its spin-down age {tsa = 
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620 kyr). PSR B0450-18 requires Vr > 360 km s'^ 
at its most probable distance inferred from bootstrap 
{D = 0.76 kpc), which may not be unreasonable, but 
in any case it has the largest uncertainties among the 
objects in our sample (c.f. !j5]). 

6.1.1. PSR B2045-16 

In our analysis of pulsar trajectories, PSR B2045— 16 
was found to pass on the sky within ~ 20 pc of the open 
cluster NGC 6604 for a radial velocity V, « -200 km s^^. 
Since the distance between the cluster and the pulsar 
was also < 20 pc, NGC 6604 could be the birth loca- 
tion of PSR B2045— 16. If it was born in this cluster, 
the kinematic age of PSR B2045-16 is Tkin ~ 1-9 Myr, 
slightly less than its spin-down age r^d — 2.85 Myr. 
NGC 6604 is an open c luster with an age ~ 6.5 Myr 
([Kharchenko et aTl I2005D . We infer that the progeni- 
tor of PSR B2045— 16 was probably more massive than 
~ 20 M0 and that the birth velocity of the pulsar was 
~ 575 km s~^. Furthermore, the mismatch between the 
spin-down age and kinematic age can be resolved if the 
initial spin period of the pulsar was Pq ~ 360 ms, assum- 
ing a fixed braking index n = 3, consistent with the popu- 
lat ion birth spin period distribution (300±150 ms) found 
by lFaucher-Giguere fc Kaspil (j2006f ) . Alternatively, if the 
initial spin period is assumed to be much smaller than 
its current period, the kinematic age implies a braking 
index n « 4, which is larger than all measured values 
to date and thus an unlikely explanation. Further analy- 
sis will need to confirm the association of PSR B2045— 16 
with NGC 6604, taking into account the proper motion of 
NGC 6604 (Ma = -0-36 mas yr"!, jis = -2.68 mas yr"! 
iKharchenko et aLll2005D . which corresponds to a motion 
on the sky of ~ 60 pc in 2 Myr. 

6.2. Electron Density Models 

Model independent distances to pulsars provide essen- 
tial calibration points for models of the Galactic electron 
densi ty distribution such as NE2001 (jCordes fc Laziol 
[2002), and one of the objectives of our observations is 
to use new pulsar distances to improve successors to 
NE2001. We have compared our astrometric distance es- 
timates with those obtained from NE20 01, as illustrated 
in Fig ure M In constructing NE2001, ICordes fc Laziol 
(|2002l ) estimated the uncertainty in the distances derived 
from the pulse dispersion measure (DM) to be 20%. For 
the purposes of our comparison, we do not include the 
NE2001 uncertainty estimate, but consider only whether 
the NE2001 distance estimate is within the 68% and 99% 
confidence intervals of the parallax distance. 

Of the 14 pulsars in our sample, the NE2001 distance 
estimates are in approximate agreement with our paral- 
lax determinations (within or just beyond the 68% con- 
fidence intervals) for 9 of the pulsars. The other five 
pulsars (B0031-07, B1508-h55, B2045-16, B1541-I-09, 
and B0450-f 55) are discussed below. 

Three of the discrepant pulsars (B0031-07, B1508+55, 
and B2045-16) have high Galactic latitudes |6| > 30° 
and low DMs DM < 20 pc cm^^. The NE2001 model 
contains a set of low electron density components de- 
signed to reproduce the effects of local Galactic structure, 
namely regions such as the Local Hot Bubble (LHB), a 
local Low Density Region (LDR), and the Local Super- 
bubble (LSB). The distances to these three pulsars are 
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Fig. 9. — Comparison of distance estimates from dispersion mea- 
sure and parallax distances for the objects in our sample as a func- 
tion of (top panel) absolute Galactic latitude \b\ and (bottom panel) 

dispersion me asure dl. Disp ersion measure distances are 

from NE2001 llCordes fc Laziol2002l) except for PSR B1541+09, for 
which NE2001 provides only a lower limit on the distance (see dis- 
cussion in ij6.2l l , and we have used an earlier distance estimate from 
ITavlor fc Cordos (1993). Error bars reflect only the uncertainties 
in the distance estimated from parallax, and PSR B0450— 18 has 
been excluded from the plot due to the large uncertainties. It 
is apparent that NE2001 underestimates distances to some of the 
pulsars in our sample at high Galactic latitudes, consistent with 
a larger scale height or irregularities in the ISM at the relevant 
scales. 

large enough that none are within any of these local, low- 
density regions. Thus, these local structures cannot be 
t he explanation f o r the systematic discrepancies. 

iGaensler et al.l ()2008 l) have re-analyzed the available 
DM data toward high Galactic latitude pulsars and find 
that the scale height of the ionized gas is approximately 
1.8 kpc, about a factor of 2 higher than the thick disk 
component in the NE2001 model. Analyzing the DM 
data in combination with emission measures (EM) de- 
rived from Ha observations, they also conclude that the 
filling factor of the ionized gas increases with distance 
above the G alactic plane. 

Using the IGaensler et al.l (|2008[ ) re-assessment for the 
ionized gas scale height, we have estimated the distances 
to these pulsars from their DMs. In general, the result- 
ing distances are slight over-estimates relative to their 
parallax distances. We conclude that the new parallax 
distances and the D Ms of these pulsars ar e consistent 
with the conclusion of IGaensler et al.l ()2008f ) of a larger 
scale height for the ionized gas, though their estimated 
scale height of 1.8 kpc may be a bit too high. A re- 
vised fit including the new distances presented here in- 
deed yields a slightly lower scale height (B. M. Gaensler, 
personal communicati on), but does n o t sign ificantly al- 
ter the conclusions of IGaensler et al] (l2008l ). We note, 
however, that the ISM is likely to be significantly irreg- 
ular or patchy on the relevant length scales, particularly 
at high Galactic latitudes, where chimmneys and voids 
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are important contributors to the structure. 

PSR B1541+09 is also at high Galactic latitude (6 = 
45°) but has a relatively large DM (35.24 pc cm'^). This 
DM is above the maximal value (« 30 pc cm~'^) from the 
thick disk component in the NE2001 model, which only 
provides lower limits on the pulsar distance. Thus, in or- 
der to reproduce the DM, the NE2001 model includes 
a "clump" of enhanced electron densi ty (contributing 
about 10% of the total DM), though iCordes fc Laziol 
([200^ could find no obvious feature along the line of 
sight that might generate this clump (e.g., an H II region 
or O star). This line of sight also passes above two spiral 
arms in the NE2001 model. One possibility is that the 
scale height of the ionized gas above one or both of these 
arms is larger than included in the model. 

Lastly, PSR B0450+55 lies at a relatively low Galactic 
latitude {b = 7?5), with DM = 14.495 pc cm^^. At the 
parallax distance to the pulsar, the DM inferred from the 
NE2001 model is 33 pc cm~'^, much larger than what is 
observed. The situation can plausibly be explained by a 
deficit of electrons in this direction. Within the NE2001 
model, there are attempts to account for such deficits 
("voids") in the electron density along specific lines of 
sight, but such modeling is limited by the sparseness of 
available DM-independent distance estimates. 

We have examined both the Virginia Tech Spectral- 
Line Survey H-g im age and the WHAM H-a survey 
(jHaffner et al.l [20031 ): neither shows any indication that 
the line of sight toward the pulsar has any deficit of 
H-a emission compared to nearby lines of sight. How- 
ever, a CO survey (Dame et al. 2001) shows a tongue of 
CO emission that crosses the line of sight to the pulsar. 
While the filling factor of the CO gas along the line of 
sight to the pulsar may not be large, its presence suggests 
that not all of the line of sight may be ionized, consistent 
with an apparent void in the electron density along the 
line of sight to this pulsar. 

6.3. Reference Frame Ties 

Pulse timing provides radio pulsar positions in the So- 
lar system reference frame, while VLBI measurements 
are tied to the distant quasars. Simply measuring pre- 
cise positions for pulsars via the two different techniques 
enables fundamental reference frame ties bet ween the So- 
lar system and the extragalactic ICRF (e.g. iBartel et al.l 
[1996). The recycled pulsar J1713-H074 7 is now one of two 
pulsars (along with PSR J0437-4715: lDeller et aL[|2008f) 
with high-precision astrometry and statistically signifi- 
cant measurements of parallax using both pulse timing 
and VLBI. In Table [4| we list the astrometric parameters 
for PSR J1713-1-0747 as determined by interferometry 
(this w ork) as well as by two independent pu lse timing 
efforts ([Splaver et al.ll2005l: iHotan et al.ll2006l) . 

The data used bv ISplaver et al.l ( 20051 ) span 12 years 
(though with a 4 year ga p), allowing compensation for 
timing noise. In contrast. iHotan et al.l {2006) used data 
spanning only 2.5 years, which did not allow corrections 
for timing noise, but allowed them to minimize systemat- 
ics through the consistency of their instrumental setup. 
We note that the overall consistency between the timing 
and VLBI results for proper motion and parallax is quite 
good, although the position differs in Declin ation at the 
~nias level. The longer time series allows iSplaver et al.l 
(|2005f ) to attain higher precision in their proper motion 



TABLE 4 
Astrometry for PSR J1713+0747 



Measurement 


VLBA 


Splaver et al. 


Hotan et al. 


Aa (s) 


0.5306(1) 


0.5307826(7) 


0.53077(1) 


AS (") 


0.519(2) 


0.52339(3) 


0.5228(2) 


fia (mas yr 


4.75i^" 


4.917(4) 


4.97(6) 


fis (mas yr-l) 


-3.67(16) 


-3.933(10) 


-3.7(1) 


TT (mas) 


0.95(6) 


0.89(8) 


1.1(1) 



Note. — Astrometric parameters for PSR J1713-I-0747 
from VLBA astrometry (this work) are compar ed to parame- 
ters derived independent ly from pulse timing iSplaver et al.l 
l2005i : [HoFan et al.|[2006D . All parameters are for equinox 
J2000; positions are shifted to match the VLBI data epoch 
of year 2002.0 (MJD 52275) and listed as offsets from 
RA 17''13™49= (Aa), Dec 07°37'47" (A<5). 

estimates, but the uncertainties in parallax are compa- 
rable between the 12 year timing data set and the 2 year 
VLBI data, illustrating the complementarity of the tech- 
niques. An ensemble of such measurements on recycled 
pulsars has great promise for reference frame ties, as well 
as for the detection of gravitational waves. 

7. CONCLUSIONS 

By simply measuring the positions of neutron stars 
over time to high accuracy, it is possible to establish con- 
straints on a variety of scientific questions. Here we have 
presented the motivation, methods, and results from a 
large VLBA astrometry program. We have described our 
initial VLA surveys to select targets and find in-bcam 
calibrators, as well as VLBA observations, pulsar gating, 
and data reduction. We have identified a few specific pit- 
falls, including proximity of the targets to the Sun, and 
the widespread presence of satellite-generated RFI. We 
have also described the use of the bootstrap method to 
comprehensively treat unmodeled systematic errors. The 
methods described here are serving as templates for our 
ongoing astrometry programs, and may serve as useful 
guidelines for the wider community. 

As a result of our VLBA astrometry program, we have 
measured new parallaxes and proper motions for 14 pul- 
sars, including tt = 0.13to;o2 mas for PSR B154H-09. 
The measurements have been exploited to investigate 
their kinematics and birth sites. Young pulsars are found 
to be leaving the Galactic plane, as expected, and we 
also find that spindown ages and kinematic ages are 
generally in reasonable agreement, with some notable 
exceptions. We have identified a plausible birth site 
for PSR B1508-I-55 in one of the Cyg OB associations 
(jChatterjee et al.l l2005| ) and for PSR B2045-16 in the 
open cluster NGC 6604 (this work). 

The new model-independent distances have also al- 
lowed us to revisit models of the Galactic electron density 
distribution. Comparing distance estimates from paral- 
lax and pulse dispersion measure, we find that NE2001 
underestimates distances for some objects at high Galac- 
tic latitudes, while for others the estimates agree with or 
are larger than the parallax distances. Our findings are 
consistent with a larger scale height or with irregularities 
in the ISM on relevant length scales. Independent dis- 
tance measurements are essential to calibrate models and 
probe structure in the ISM, and the present sample will 
be incorporated into future iterations of the electron den- 
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sity distribution model, ultimately improving DM-based 
distance estimates for all pulsars. 

Finally, we have compared our precise astrometry on 
the recycled pulsar J1713+0747 with two independent 
results from pulse timing. The comparison can be used 
to verify consistency between the extragalactic ICRF and 
the Solar system reference frame, especially as part of an 
ensemble of recycled pulsars. 

The Very Long Baseline Array is the only instrument 
that offers full-time, dedicated VLBI capabilities. Our 
results once again illustrate the power of the instrument 
to provide astrometry even at lower frequencies with high 
precision and excellent accuracy. 
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